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Abstract. We present high-resolution two-dimensional velocity fields 
in Ha and CO of the nearby dwarf galaxy NGC 2976. Our observations 
were made at both higher spatial resolution (~ 75 pc) and higher velocity 
resolution (13 km s _1 in Ha and 2 km s _1 in CO) than most previous 
studies. We show that NGC 2976 has a very shallow dark matter density 
profile, with p(r) lying between p oc r~ ' 3 and p oc r . We carefully test 
the effects of systematic uncertainties on our results, and demonstrate 
that well-resolved, two-dimensional velocity data can eliminate many of 
the systematic problems that beset longslit observations. We also present 
a preliminary analysis of the velocity field of NGC 5963, which appears 
to have a nearly NFW density profile. 



1. Introduction 

It is well-known by now that there is a substantial disagreement between the 
observed dark matter density profiles of many dwarf and low-surface bright- 
ness galaxies and the density profiles predicted by numerical Cold Dark Matter 
(CDM) simulations (e.g., Flores & Primack 1994; Burkert 1995; Navarro, Frenk, 
k White 1996, hereafter NFW; Moore et al. 1999). The significance of this 
disagreement, though, remains controversial. A number of authors attribute the 
problem to failures of the simulations, or of the CDM model itself (de Blok et al. 
2001a; de Blok, McGaugh, & Rubin 2001b; Borriello & Salucci 2001; de Blok, 
Bosnia, &: McGaugh 2003), while others argue that systematic uncertainties in 
the observations make such conclusions premature (van den Bosch et al. 2000; 
van den Bosch & Swaters 2001; Swaters et al. 2003, hereafter SMVB). 

We address this controversy with a new study that combines a number of 
techniques to overcome the systematics in the observations. We are acquiring 
very high-quality data on a limited sample of galaxies to investigate the impor- 
tance of systematic effects in detail. The results of this study should make clear 
whether systematic problems in the data are at fault, or whether there actually 
is a fundamental conflict between the theory and the observations. 

Our program includes (1) two-dimensional velocity fields obtained at opti- 
cal (Ha), millimeter (CO), and centimeter (H i) wavelengths, (2) high angular 
resolution (~ 5"), (3) high spectral resolution (<; 10 km s _1 ), (4) multicolor 
optical and near-infrared photometry, and (5) nearby dwarf galaxies as targets. 
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The combination of these features greatly reduces our vulnerability to systematic 
uncertainties (see Simon et al. 2003). 

In two previous papers, we reported on rotation curve studies of the dwarf 
spiral galaxies NGC 4605 and NGC 2976 (Bolatto et al. 2002; Simon et al. 2003). 
In this paper we highlight some of those results, focusing on tests for systematic 
uncertainties in the NGC 2976 data set, and we also present preliminary results 
on the third galaxy in our study, NGC 5963. 

2. The Dark Matter Halo of NGC 2976 

NGC 2976 is a regular Sc dwarf galaxy located in the M81 group at a distance 
of 3.45 Mpc. The galaxy has absolute magnitudes of Mg = —17.0 and Mx = 
—20.2, and a total mass of 3.5 x 1O 9 M ; it is somewhat less luminous and less 
massive than the Large Magellanic Cloud. In optical and near-infrared images it 
is clear that NGC 2976 is a bulgeless, unbarred, pure disk system, which makes 
it an ideal galaxy for mass modeling. 

2.1. Observations 

Our observations of NGC 2976 include a two-dimensional Ha velocity field 
(obtained with a multi-fiber spectrograph on the WIYN telescope), a two- 
dimensional CO velocity field (obtained with the BIMA interferometer), multi- 
color optical photometry, and near-infrared 2MASS imaging. The velocity fields 
are both shown in Figure 1. For details of the data reduction and analysis, see 
Simon et al. (2003). 

2.2. Rotation Curve 

We used various tilted-ring modeling algorithms to convert the observed veloc- 
ity field into a rotation curve (see Simon et al. 2003). The rotation curve of 
NGC 2976 is well-described by a single power law from the center of the galaxy 
out to a radius of almost 2 kpc, as displayed in Figure 2a. The rotation curve 
only begins to deviate systematically from power law behavior at r ~ 110" (1.84 
kpc). The total (baryonic plus dark matter) density profile corresponding to the 
rotation curve is Ptot = 1.6(r/l pc) _0 ' 27±0 ' 09 Mq pc~ 3 . Because the baryons 
are almost certainly centrally concentrated, this density profile represents the 
cuspiest possible shape for the dark matter halo. 

2.3. Dark Matter Density Profile 

To explore the full range of possible dark matter density profile slopes, we must 
account for the rotational velocities contributed by the baryonic components of 
the galaxy. Because our images of NGC 2976 do not reveal a bulge or a bar, and 
its nucleus is dynamically unimportant, the only relevant reservoirs of baryons 
to consider are the stellar and gaseous disks. We calculated the rotation curves 
due to the stars and gas directly from the observed surface density profiles, 
assuming an infinitely thin disk. The only free parameter in these calculations 
is the K-band mass-to-light ratio of the stellar population (M*/Lk)- 

Under the assumption that the density profile can be described with a power 
law, p DM oc r~ a ° M , we found that if M^/Lk > 0.19M Q /L Q x, «dm < and the 
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Figure 1. (a) Ha velocity field of NGC 2976 from WIYN observa- 
tions. The contours represent integrated Ha intensity, (b) CO velocity 
field from BIMA observations. The contours represent integrated CO 
intensity. The angular resolution of each dataset is shown in the upper 
left corners. 



density of the dark matter halo is increasing with radius. Such a dark matter con- 
figuration is probably unphysical, so we consider 0.19 Mq/Lqk to be a firm up- 
per limit to the stellar disk mass-to-light ratio, with the corresponding lower limit 
to a DM of 0. Because of the extremely low value of the maximal disk mass-to- 
light ratio, the galaxy must contain an essentially maximal disk. The dark mat- 
ter density profile for the maximal disk is p DM = 0.1(r/l pc)~ o,ol=l:O ' 13 M0 pc~ 3 . 
This represents the most likely shape for the dark matter halo. Since the slope 
of the total density profile of the galaxy represents the absolute upper limit for 
the slope of the dark matter density profile, we conclude that the dark matter 
density profile is bracketed by p DM oc r -°- 27±0 - 09 an d p DM oc r°. 



3. Systematics 

Several authors have recently discussed the systematic uncertainties that can 
significantly alter observed rotation curves. SMVB used simulated observations 
to argue that systematic effects alone could account for the difference between 
the predicted and observed density profile slopes, de Blok et al. (2003) car- 
ried out similar simulations, but concluded that the systematic effects were not 
strong enough to cause cuspy density profiles to appear to have cores. We have 
shown that, for the specific case of NGC 2976, the known systematic problems 
do not distort the rotation curve (Simon et al. 2003; this work). The worst 
systematics can probably be minimized or avoided in general by using two- 
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Figure 2. (a) Minimum disk rotation curve of NGC 2976. The ob- 
served rotation velocities are plotted as black circles and the error bars 
are combined statistical and systematic uncertainties. The dashed and 
dotted curves represent the rotation velocities due to H I and H2, re- 
spectively. A power law fit to the rotation curve is shown by the solid 
black curve. The corresponding density profile is p oc r~°' 27 . Residuals 
from the fit are displayed in the lower panel, and la and 2a depar- 
tures from the fit are represented by the shaded regions, (b) Maximum 
disk rotation curve of NGC 2976. The stellar disk (thick gray curve) 
is scaled up as high as the observed rotation velocities (gray circles) 
allow. The stellar disk shown here has M^/Lk = 0.19 Mq/Lqk- Af- 
ter subtracting the rotation velocities due to the stars and the atomic 
and molecular gas from the observed rotation curve, the dark mat- 
ter rotation velocities are displayed as black circles. The two missing 
data points near the center of the galaxy had v ro t < v*, r otj yielding 
imaginary Vhalo- The solid black line is a power law fit to the halo ve- 
locities (for 14" < r < 105") which corresponds to a density profile of 
Pdm of r~ - 01 . The halo residuals after the power law fit are displayed 
in the bottom panel. 



dimensional velocity fields and by making velocity measurements at very high 
precision (Blais-Ouellette et al. 1999; van den Bosch & Swaters 2001; Bolatto 
et al. 2002; SMVB), but additional observations confirming this expectation for 
more galaxies are needed. 

3.1. Rotation Curve Fitting Systematics 

Because of the high precision of our velocity measurements, the statistical un- 
certainties on both the rotation curve and the radial velocity curve are negligible 
(less than 1 km s^ 1 everywhere). Therefore, the limiting factors on the accuracy 
of the rotation curve are the systematic uncertainties associated with our fit. 
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Uncertainties in Geometric Parameters Before computing the uncertainties 
on the rotation velocities themselves, we must first determine the uncertainties 
on each of the parameters that are used to calculate the rotation velocities: 
the center, position angle (PA), inclination, and systemic velocity. We used a 
bootstrap resampling technique to calculate the value and uncertainty of these 
parameters. We constructed 200 bootstrap samples of the velocity field and 
ran our tilted-ring modeling routine (ringfit) on each of them to determine 
the best-fit kinematic center. We then defined the center of the galaxy to be 
the median of these 200 measurements, and the systematic uncertainty on the 
center position to be the dispersion of the measurements. We found that the 
location of the kinematic center of NGC 2976 is consistent with the optical 
nucleus, with an uncertainty of 2" in both a and 5. We used the same bootstrap 
method to measure the kinematic PA of the galaxy and its uncertainty, finding 
that the kinematic PA matches the photometric one and has an uncertainty of 
5°. It is not possible to determine a kinematic inclination angle for NGC 2976 
because the rotation curve is too close to solid-body. However, the photometric 
inclination angle is well determined, and we estimate that the uncertainty is 3°. 

The uncertainty in the systemic velocity was calculated in a different way. 
The systemic velocity of each ring was left free to vary during our tilted-ring 
fitting, resulting in 23 independent measurements. The overall uncertainty from 
these measurements is 1.8 km s -1 . Since v sys is always a free parameter, this 
uncertainty does not directly factor in to the uncertainty on the rotation curve. 

Uncertainties in Rotation Velocities and Radial Velocities Using the measured 
uncertainties in the center position, PA, and inclination angle, we calculated the 
resulting uncertainties on the rotation velocities and the radial velocities with a 
Monte Carlo study. We generated 1000 random centers, PAs, and inclinations, 
assuming a Gaussian distribution for each of the parameters, and ran ringfit 
with each set of parameters. The standard deviation of the 1000 rotation veloc- 
ities in each ring was defined to be the systematic error of that rotation velocity, 
and the systematic errors in the radial velocities and systemic velocities were 
calculated in the same way. The systematic errors on the rotation curve range 
from 2.1 km s _1 to 5.5 km s _1 . The sum in quadrature of these errors and the 
statistical errors translates to a total uncertainty on the density profile slope a 
of 0.09. We conclude that, even in the minimum disk case, an NFW or steeper 
density profile in NGC 2976 is strongly ruled out. 

3.2. Other Systematic Effects 

In some studies, beam-smearing is one of the dominant systematic effects. The 
velocity field of NGC 2976, though, is extremely well-resolved — our spatial 
resolution is less than 100 pc in both CO and Ha. The measured velocity 
gradient is small (~ 50 km s _1 kpc -1 ) and the constant-density core that we 
detected is ~ 50 resolution elements across, indicating that beam-smearing is not 
responsible for this result. We also calculated the asymmetric drift correction 
to the rotation curve. Its effect is to increase the maximum disk mass-to-light 
ratio, and to make the rotation curve slightly more linear, but the density profile 
does not change significantly. 



6 



Simon et al. 



3.3. Comparing Velocities Derived Prom Different Tracers 

Some recent studies in the literature have shown that, beam smearing questions 
aside, there do not appear to be systematic offsets between H I and Ha rotation 
velocities (e.g., McGaugh, Rubin, & de Blok 2001; Marchesini et al. 2002). With 
a handful of exceptions, though, these studies employed longslit Ha data, so the 
comparisons essentially took place only along the major axis. In addition, the 
spatial and velocity resolution of the H I and Ha data were often quite different. 

Simon et al. (2003) presented for the first time the data necessary for a two- 
dimensional comparison across a dwarf galaxy of the CO and Ha velocity fields. 
The angular resolution of the two datasets is similar (6" and 4", respectively), 
and although the CO velocity resolution is better by a factor of ~ 6, the higher 
signal-to-noise at Ha enables us to measure the velocities with comparable pre- 
cision. We constructed a unique one-to-one mapping between the two velocity 
fields. The mean difference between v Ha and v co is 1 km s~ x and the rms is 
5.3 km s -1 , with the comparison being made at 173 independent points. Sim- 
ilar studies in the Milky Way found that the dispersion between the velocities 
of molecular clouds and the associated Ha-emitting gas was 4-6 km s -1 (Fich, 
Dahl, & Treffers 1990), so much of the scatter we observe between the two in 
NGC 2976 may be intrinsic to the process of H n region formation rather than 
caused by observational uncertainties. 

This comparison strongly suggests that both the CO and the Ha data accu- 
rately reflect the gravitational potential of the galaxy, and neither is significantly 
compromised by systematic effects. Extinction, for example, is not distorting the 
Ha velocity field. If this result can be shown to hold for a few additional galaxies, 
then it will no longer be necessary to observe more than one velocity field tracer 
per galaxy. One further test that should be carried out is to compare the velocity 
field from a gaseous emission-line tracer to that from a stellar absorption-line 
tracer to assess the agreement between the stellar and gas kinematics. Such stud- 
ies have been done before (e.g., Mulder 1995; Bottema 1999), but usually only 
with long-slit observations. Obtaining two-dimensional absorption-line velocity 
fields is currently observationally challenging but feasible. 

3.4. Noncircular Motions 

Perhaps the most surprising result from Simon et al. (2003) was the finding that 
NGC 2976 contains strong radial motions. In optical and near-infrared images, 
the galaxy is very regular, with no hint of a bar or any other deviation from 
axisymmetry that might affect the kinematics. The finding that the velocity 
field is clearly different from a purely rotating disk was therefore unexpected. 
Near the center of the galaxy, the direction of the maximum velocity gradient is 
up to ~ 40° away from the photometric major axis. 

The velocity data can be adequately described by a model in which the 
galaxy is undergoing circular rotation, but only if the kinematic PA declines 
monotonically from ~ 6° to ~ —37° over the central kiloparsec. As mentioned 
above, though, this position angle variation is inconsistent with the photometry. 
The alternative — and more likely — model is one in which NGC 2976 contains 
substantial radial flows. The origin of such motions, however, is unclear. 

There are several plausible mechanisms for creating radial motions. An 
intriguing possibility is that the radial motions could be a result of the dark 
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matter halo having a triaxial rather than spherical shape. CDM halos are ex- 
pected to be moderately triaxial (e.g., Dubinski & Carlberg 1991; Warren et 
al. 1992), and the velocity field of a galaxy embedded in a triaxial halo would 
exhibit noncircular motions. However, since the details of such a velocity field 
have not yet been investigated with simulations, we cannot compare our results 
to theoretical predictions. We plan to carry out simulations of the kinematics of 
a disk within a triaxial halo to see how the halo shape influences the gas kine- 
matics and whether the effect is consistent with the velocity field of NGC 2976. 
This could become a new technique for measuring the triaxiality of galaxy dark 
matter halos. 



4. The Dark Matter Halo of NGC 4605 

NGC 4605 is another Sc dwarf galaxy, located at a distance of 4.3 Mpc. Like 
NGC 2976, it appears to be a pure disk system. We used a long-slit Ha spec- 
trum and a two-dimensional CO velocity field to measure its rotation curve 
(Bolatto et al. 2002). We found that the galaxy probably contains a maximal 
disk, and that its dark matter density profile goes as p oc r~ 65 (see Fig. 3), 
intermediate between a core and a cusp. The quality of our original observations 
of NGC 4605 precluded an investigation of the systematics, but we are in the 
process of analyzing new observations to study their impact. 

5. The Dark Matter Halo of NGC 5963 

NGC 5963 is a more massive (~ 1O 1O M0) and more distant (~ 10 Mpc) galaxy 
than NGC 2976 and NGC 4605. Its surface brightness profiles are more com- 
plicated than those galaxies as well, showing bright emission at the center, sur- 
rounded by a transition region in which the surface brightness falls off extremely 
rapidly (scale length ss 300 pc), and then a very low surface brightness expo- 
nential disk extending from 1.5 kpc out to beyond 3 kpc. The interpretation of 
this behavior near the center of the galaxy is not straightforward; the galaxy is 
definitely not strongly barred, but we cannot rule out a bulge component. 

5.1. Observations 

We obtained a two-dimensional Ha velocity field of NGC 5963 with the WIYN 
telescope, using the same technique as for NGC 2976. The spatial resolution 
of these data is 190 pc and the velocity resolution is 13 km s _1 . We detected 
Ha emission out to a radius of ~ 3 kpc. We also obtained a two-dimensional 
CO velocity field from BIMA. The CO emission in this galaxy is limited to the 
central kiloparsec, but adding the CO data improves our ability to constrain the 
innermost part of the density profile. 

5.2. Rotation Curve and Stellar Disk 

The combined CO and Ha rotation curve of NGC 5963 is displayed in Figure 4. 
It is immediately clear that the density profile of NGC 5963 has a very different 
structure from those of NGC 2976 and NGC 4605. Instead of a mostly linear 
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Figure 3. (a) Minimum disk rotation curve of NGC 4605. The black 
curve is a p oc r -1 ' 1 power law fit, which agrees quite closely with 
the rotation curve beyond 25". At smaller radii, though, it signifi- 
cantly overestimates the rotation velocities. A separate inner power 
law (p oc r~ ' 4 ) is necessary to fit the entire rotation curve. The bot- 
tom panel shows the residuals from the fit, with la and 2a departures 
represented by the shaded regions. Statistical error bars are displayed 
in the residual plot, (b) Maximum disk rotation curve of NGC 4605. 
The light gray circles represent the observed rotation curve, and the 
black circles represent the rotation curve of the dark matter halo after 
the stellar disk contribution (thick gray curve) has been subtracted. 
The black curve is a p oc r~°- 65 power law that provides a good fit to 
the dark matter rotation curve at all radii. 



increase with radius, the rotation curve of NGC 5963 has the classic spiral galaxy 
shape with a steep inner rise and a nearly flat outer portion. 

We perform isophotal fits to the I-band image of NGC 5963 to obtain the 
shape of its surface brightness profiles. We then calculate the rotation curve due 
to the stellar disk numerically, as described before. The resulting stellar disk 
rotation curve can be seen in Figure 4. The maximum disk limit for this galaxy 
is set by the innermost point of the rotation curve. We derive a maximum disk 
mass-to- light ratio of M^/Lj = 0.7 'Mq/L & j. Even the maximum disk is unable 
to account for the observed rotation speed of the galaxy beyond a radius of ~ 6". 
The calculated value of the maximum disk mass-to-light ratio is not unusually 
low, so a maximum disk is not required in NGC 5963. The disk can be less 
massive than shown in Figure 4, but since the galaxy is dark matter-dominated 
the structure of the dark halo does not depend strongly on the assumed stellar 
mass-to-light ratio. 
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Figure 4. (a) Minimum disk rotation curve of NGC 5963. The black 
curve shows an NFW fit to the rotation curve with a concentration 
of 20. Residuals from the fit are plotted in the bottom panel. Un- 
like the previous two galaxies, a power law is not a very good fit and 
the NFW form fits much more closely, (b) Maximum disk rotation 
curve of NGC 5963. Even after removing a maximal stellar disk (thick 
gray curve; M^/Lj = O.7M /L /), the rotation curve due to the dark 
matter halo is well-fit by an NFW potential. 



5.3. Dark Matter Density Profile 

Fitting a single power law to the rotation curve of NGC 5963 does not produce 
a very good fit. The fit is substantially improved by using separate inner and 
outer power laws, corresponding to an inner density profile of p oc r -11 and a 
steeper outer density profile of p cc r~ 1A . An even better fit can be obtained by 
fitting an NFW rotation curve to the data; the reduced % 2 value of this fit is only 
1.1. Although the concentration parameter of the halo is not well-constrained, 
the scale radius is about 3 kpc and the halo has V200 ~ 90 km s -1 . An NFW 
fit can accurately describe the rotation curve for any stellar mass-to-light ratio 
(see Fig. 4) . This galaxy may be the first low- mass system found that contains 
a density profile that matches the predictions of the numerical simulations. 



6. Conclusions 

We have used high-resolution two-dimensional velocity fields to study the dark 
matter density profiles of NGC 2976, NGC 4605, and NGC 5963. We showed 
that these three galaxies contain very different density profiles. For NGC 2976, 
we presented one of the most detailed velocity fields of a dwarf galaxy outside 
the Local Group. Even with these high-quality data, we found that NGC 2976 
contains a constant-density core and an NFW halo is strongly ruled out. The 
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dark matter halo of NGC 4605 is intermediate between a constant-density core 
and a cusp, and NGC 5963 has a cuspy halo that is almost perfectly consistent 
with an NFW profile. 

We conclude that some galaxies do not contain central cusps, even when 
systematic uncertainties are accounted for as carefully as possible. We also point 
out that we find no evidence to support the prediction that all dark matter 
halos share a universal shape. If these results hold as our sample grows, it will 
demonstrate that the conflict between the observations and the simulations is 
not caused by systematic problems with the observations. Instead, more effort 
may be needed to investigate what could be missing from the simulations that 
would cause them to overestimate density profile slopes. 
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